#ifndef GAUSS_HPP_
#define GAUSS_HPP_

/*! \file gauss.hpp
 * \brief declaration of Gauss points as global variables
 */

extern double xi[51][51]; //!< Legendre Gauss points 
extern double wi[51][51]; //!< Legendre Gauss weights 
extern double Pi; //!< Pi

/*! \fn void init_gauss ()
  \brief has to be called once to initialize the global variables of Gauss points
 */

void initGauss(double *wi, double *xi, int _GaussDegree);



#endif // GAUSS_HPP_ 
